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ABSTRACT 



Aims. A major goal of our survey is to significantly increase the number of detected extra-solar planets in a magnitude-limited sample 
to improve our knowledge of their orbital elements distributions and thus obtain better constraints for planet-formation models. 
Methods. Radial- velocity data were taken at Haute-Provence Observatory (OHP, France) with the ELODIE echelle spectrograph. 
Results. We report the presence of a planet orbiting HD 196885 A, with an orbital period of 1349 days. This star was previously sug- 
gested to host a 386-day planet, but we cannot confirm its existence. We also detect the presence of a stellar companion, HD 196885 B, 
and give some constraints on its orbit. 

Key words, stars: individual: HD 196885 A - stars: individual: HD 196885 B - stars: binaries: visual - stars: planetary systems - 
techniques: radial velocities - methods: observational 



1. Introduction 

The ELODIE Planet Search Survey was an extensive radial- 
velocity northern survey of dwarf stars at the Haute-Provence 
Observatory (OHP, France) u sing the ELODIE hi gh-precision 
fiber-fed echelle spectrograph (Baran ne et al.lll996h mounted on 
the Cassegrain focus of the 1 .93-m telescope. It started to operate 
at the end of 1993 and acquired data until the summer of 2006, 
when it was replaced by the new e chelle spectrograph SOPHIE 
dBouchv & The Sophie Teamll2006l) . 

This survey is part of a large effort aiming at an extra-solar 
planet search through radial-velocity measurements, in order to 
characterize the zoo of exoplanets and to bring strong constraints 
on their processes of formation and evolution. ELODIE was re- 
sponsible for the findings of man y extra-solar planets, a mong 
them the first hot Jupiter, 5 1 Peg b (iMavor & Ouelozll 19951) . The 
original sample consisted of 142 stars, but a new sample of 330 
stars was defined in 1997. Details ab out the progr am and the 
surveyed sample can be found in lPerrier et al.l ([2003). 

In this paper we present the detection of two bodies around 
the star HD 196885 A, one in the mass range of planets and the 
other believed to be a stellar companion. This star was previously 
described to harbor a planet with a period of about 386 days (al- 
though no published reference is known), but we could not con- 
firm its presence in our observations (our analysis shows a planet 
with an orbital period of about 1349 days instead). The stellar 
companio n was also recently observe d using NACO adaptative 
optics by IChauvin et al.l (2006, 2007). The stellar properties of 
HD 196885 A are briefly recalled in Sect. 2 and the radial veloc- 
ities with the two detected companions are described in Sect. 3. 



In Sect. 4 we discuss the possibility of existence of more com- 
panions, in particular around 368 days. Last section is devoted 
to our conclusions on this system. 

2. HD 196885 A stellar characteristics 

HD 196885 A was observed by the HIPPARCOS astrometric 
satellite (HIP 101966). A high-pre cision spectroscopi c study of 
this star was also performed by ISousa et at] (|2006) in order 
to examine the metallicity distribution of stars hosting planets. 
Moreover, this star was also studied by the near-infrared survey 
with adaptive optics of faint circumstellar environments, sensi- 
tive to companions w i thin th e stellar and the sub-stellar domains 
dChauvin et al.l [20061 120071) . Observed and inferred stellar pa- 
rameters from these different sources are summarized in TableQ] 
In the HIPPARCOS catalogue, HD 196885 A is given a spec- 
tral type F8 IV, a visual magnitude V = 6.398 and a color in- 
dex B - V = 0.559. The measured parallax (30.31 ± 0.81 mas) 
leads to a distance of 33.0 + 0.9 pc and an absolute magnitude 
of My = 3.8 ± 0.1. From CORALIE spectra. ISousa et alj (2006) 
derived an effective temperature T e n = 6340 ±39 K, a gravity log 
g = 4.46 ± 0.02 and a high metal content [Fe/H] = 0.29 + 0.05 
(TableQ}. The bolom etric correction (BC = -0.006) is computed 
from Flower (1996) using the spectroscopic T e ^ determination. 
The bolometric magnitude is then Mgoi = 3.791, which allow 
us to derive a stellar luminosity of L — 2.40 L Q . ISousa et all 
( 20061) also compute the mas s, M — 1.33 M B , from evolutionary 
tracks using Geneva models CSc haller et al.ll 99"2t ISchaerer et al.l 
119931) . From the B - V value and ELODIE correlation functions 
we find vsini = 7.3 + 1.5 [km/s], meaning that the star is ro- 
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Table 1. Observed and inferred stellar parameters for 
HD 196885 A. Photometric and spectral type are from 
Astrometric parameters are from 



HD 196885A 



ELODIE + CORALIE 



IChauvin et al. 


(I2006I). 


HIPPARCOS ( 


ESaT|1997) 



parameters T e ^, log g and [Fe/H] are from ISousa et al. (2006). 
The bolometric correction is computed from iFlowerl (l996) 
using the spectroscopic T e ff determination . The a ctivity level and 
the rotation period are from lWright et alj d2004l) . The given age 
was o btained following the Bayesian approach of iPont & Ever! 
(200%. 



Parameter 




HD 196885 A 


Spectral Type 




F8V 


V 




6.398 


B-V 




0.559 ± 0.006 


n 


[mas] 


30.31 ±0.81 


d 


[pc] 


33.0 ± 0.9 


M v 




3.8 + 0.1 


BC 




-0.006 


M Bo i 




3.79 


L 


[i Q ] 


2.40 


[Fe/H] 




0.29 ± 0.05 


log*™ 




-5.01 




[days] 


15. 


M 


[M ] 


1.33 


Teg 


[K] 


6340 ± 39 


l°gg 


[cgs] 


4.46 ± 0.02 


v sin i 


[km/s] 


7.3 ± 1.5 


age 


[Gyr] 


2.0 ±0.5 



tating fast. Howeve r, its measured act ivity level is very weak: 
\ogR' HK = -5.01 (IWright et al.l 12004). From this value, these 
au thors inferred a rotatio n period of about 15 days. According 
to lPace & Pasquinil (120041) . the activity level becomes weak and 
constant after ~1.5 Gyr, setting an inferior limit on the age of 
the star. O n the other hand, fol lowing the Bayesian approach de- 
scribed in Pont & Ever (2004) we estimate a maximum age of 
2.5 Gy. The derived stellar atmospheric parameters and age are 
compatible with the expected values for a high metallicity late- 
F dwarf (F8 V), and are slightly at odds with the evolutionary 
status given by the HIPPARCOS catalogue. 

3. Orbital solutions for the HD 196885 system 

The ELODIE observations of HD 196885 A started in June 1997 
and the last data acquired are from August 2006, since the 
ELODIE program was closed shortly after that date. The peculiar 
variations of the radial velocities (FigUJ and also non-confirmed 
announcements from other research teams (see section |4)l, pre- 
vented us from announcing this system earlier. However, the re- 
cent detec tion of a visual small stella r companion close to the 
main star dChauvin et al.ll2006l 120071) . confirmed our suspicion 
of a long term drift of the radial velocities. Superimposed with 
the drift we can also observe a regular variation of a few years 
signalizing the presence of a sub-stellar companion. 

Before the ELODIE program, the star HD 196885 A had al- 
ready been followed betw een June 1982 a nd August 1997 by the 
CORAVEL spectrometers (Baran ne et al. 1 119791) mounted on the 
1-m Swiss telescope at Haute-Provence Observatory and on the 
1.54-m Danish telescope at La Silla Observatory (ESO, Chile). 
The precision of CORAVEL is ~0.3 Km/s, not enough to detect 
planetary objects, but very useful to help us to constrain the orbit 
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Fig. 1. ELODIE and CORALIE radial velocities for HD 196885 A, 
superimposed on a two-keplerian orbital solution (Tabf2|i. 



of the stellar companion. From April 1999 to November 2002 a 
series of radial velocity measur ements were also tak en using the 
CORALIE echelle spectrograph dOueloz et alJl2000l) mounted on 
the 1.2-m Swiss telescope at La Silla. This four year observa- 
tional sequence is important to confirm the presence of the plan- 
etary companion, since the precision of CORALIE is slightly bet- 
ter than the precision of ELODIE. 

With 111 radial-velocity measurements (69 from ELODIE, 9 
from CORAVEL and 33 from CORALIE), spanning -14 years of 
observations, we are able to describe the orbit of the sub-stellar 
body in the system, as well as slightly constrain the orbit of 
the stell ar companion. U sing the iterative Levenberg-Marquardt 
method dPress et alJll992l) . we first attempt to fit the complete 
set of radial velocities with a single orbiting companion and 
a quadratic drift. This fit yields a planetary companion with 
P = 1347 days, e = 0.44, a minimum mass of 2.9Mj up and 

an adjustment of y[x* = 1.811 and rms = 14.61 m/s. We then 
fit the radial velocities using a model with two Keplerian orbits 
(FigsQ]and[2]). It yields for inner planet P = 1349 days, e = 0.46 
and a minimum mass of 3.0Mj up , while for the outer compan- 
ion we have P ^ 20000 days, e = 0.41 and a minimum mass 
of 0.34 M Q (Tabf2]i. Despite all the uncertainties in the orbital 
parameters, the use of a Keplerian orbit for the massive outer 
body in the system proved to be a good approach, better than the 
quadratic drift, since the reduced -^)^ is now 1 .494 and the ve- 
locity residuals drop to rms = 1 1 .87 m/s (Fig|2|i. Of course one 
can always argue that increasing by four the number of free pa- 
rameters in the model will improve the fit, but at least we are 
able to give constraints for the orbital parameters of the outer 
body. For instance, we find that it must present an orbital pe- 
riod P > 40 yr, a semi-major axis a > 14 AU and a minimum 
mass M sin i > 0.28 M , since all other fits coherent with the 
data found larger values for these parameters (Tabf3]l. 

As expected, the orbital parameters of the outer body still 
present some uncertainties around the best fitted value. This is 
particularly true for the orbital period, meaning that this pa- 
rameter may assume rather different values. The Levenberg- 
Marquardt method converges to a minimum^ 2 , but other close 
local minima may represent as well a good fit to our data. For 
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Table 2. Orbital parameters for two bodies orbiting 
HD 196885 A, obtained with a two-keplerian fit to observa- 
tional data. Errors are given by the standard deviation cr. The 
orbital period of the outer body is much longer than the data 
acquired so far (-14 yr) and thus we are unable to completely 
constrain its orbit. However, we noticed that it is actually better 
to fit a complete elliptical orbit to the present data, rather than 
use a quadratic drift, for which we obtained y/x* — 1.811 and 
rms = 14.61 m/s. Alternative solutions for the outer body with 
longer orbital periods exist, that match the data almost as good 
as the best fit (Tab©. 



Param. 


[unit] 


HD 196885 Ab 


HD 196885 B 


V 


[km/s] 


-31.666 


+ 0.455 


P 


[year] 


3.69 ± 0.03 


55.43 ± 19.48 


e 




0.462 ± 0.026 


0.409 ± 0.038 


OJ 


[deg] 


91.4 ±4.1 


227.6 ± 23.4 


K 


[m/s] 


40.48 ± 2.32 


2041 ± 270 


T 


[JD-2400000] 


51236 + 18 


45928 ± 2638 


a\ sin ( 


[1(T 3 AU] 


4.45 


3466 


f(M) 


[10- 3 M Q ] 


6.5 x 10- 6 


13.56 


M sin i 


[M Jup ] 


2.96 


351.5 


M sin i 


[M Q ] 


0.0028 


0.34 
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2.63 


17.23 
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[m/s] 


11.87 


V? 
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Fig. 2. Phase-folded radial velocities measurements and best fit 
for the planetary companion of HD 196885 A (top) and the its 
stellar companion (bottom). For each body the contribution by 
the other companion has been subtracted from the observational 
data. The orbital period of the inner body is P = 3.7 yr, while 
for the outer body we have P — 55 yr (Tab©. However, the fact 
that we are unable to completely cover in phase the orbit of the 
stellar companion, clearly indicates that its orbital period may be 
considerably longer (Tab©. 



a particular choice of the orbital period, all the other parameters 
can adjust in a way such that the x 2 value remains close to the 
value from the best fit (Tab©. Thus, any orbital solution with a 
X 2 and an rms close to the values from Table [2] can be seen as a 
good model of the system (eg. Tab©. 



Table 3. Alternative orbital parameters for the stellar companion 
HD 196885 B. With the present ~14 years of observational data 
it is still impossible to co mpletely constrain th e outer body or- 
bit. Observations done by Chauvi net alj (|2007) point to a small 
stellar companion with 0.5 - 0.6 M at about 23 AU. Fixing the 
orbital period at larger values, we are able to find orbital solu- 
tions compatible with those observations. 
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a 
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OJ 


T 
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[yr] 


[AU] 




[deg] 


[JD] 


[M ] 


[m/s] 




40 


13.7 


0.50 


204 


45983 


0.28 


11.95 


1.504 


60 


18.3 


0.41 


233 


45916 


0.36 


11.87 


1.494 


80 


22.4 


0.42 


250 


46136 


0.43 


11.88 


1.495 


100 


26.6 


0.45 


258 


45987 


0.55 


11.88 


1.495 


120 


30.4 


0.49 


265 


46110 


0.61 


11.89 


1.496 



IChauvin et ail (120061 l2007h determined that HD 196885 is a 
close visual binary. They derive a possible mass for the compan- 
ion of 0.5 - 0.6 M Q at a projected physical distance of -23 AU 
and spectral type Ml V. As they noticed, for those values of the 
mass and semi-major axis we should be able to monitor its or- 
bital motion using radial velocity measurements of the main star. 
In order to test if our drift can be interpreted as the observed 
stellar companion, we fitted our data by fixing orbit al periods 
compared to those forecasted by IChauvin et alj (120071) . As it is 
shown in Table [3] we are able to find multiple possibilities that 
match their observations. Notice also that the values of the mass 
listed in Table [3] correspond to minimum masses. Because of 
the factor sin i, the real masses can be larger. We then conclude 
that the outer companion of the star HD 19 6885 A is almost cer- 
tainly the stellar mass object observed by Cha uvin et al.l ([2006, 
2007), henceforward denoted by HD 196885 B. The inner com- 
panion of HD 196885 A remains a planetary body, henceforward 
HD 196885 Ab, with orbital parameters given by those listed in 
Table© 

Finally, we note that the HD 196885 system may be a 
part of a wider binary with t he star BD+ 104251 B located 
-191.9" North of HD 19688 5 A dTrilling et alj|2007l) . At a dis- 
tance of 33 pc measured bv lPerrvman et a 1.1 (119971) . we estimate 
a semi-major axis of 6330 AU, that is, a drift contribution of 
-1 cms^'yr -1 to the radial velocity of HD 196885 A (assuming a 
circular orbit), and thus impossible to be detected with ELODIE. 



4. A second planet around one year? 

In 2004 it was reported a P = 386 days companion to the 
HD 196885 A star, by the California & Carnegie Planet Search 
Team on their webpageQ detected through radial-velocity mea- 
surements obtained with the Lick survey. T his planet is also 
listed i n the the exoplanets E ncyc lopedia^ , by iFischer & ValentH 
(120051) . IChauvin et all d2007l) and lMarchil(l2007l) . Unfortunately 
it is not known any publicati on describing the system as no- 
ticed by many other author s (Chauvin et al. 2006; Sousa et al.l 
2006; iBonavita & Desideral 120071; iDesidera & Barbieril [2007b . 
Moreov er, the system is no t reported in the 2006 exoplanet cat- 
alog by iButler etldl (120061) . so we assume there was no confir- 
mation of this planet after the first announcement. 



http://exoplanets.org/esp/hdl96885/hdl96885.shtml 
http://exoplanet.eu/planet.php ?pl=HD+196885&p2= 
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Fig. 3. Frequency analysis (top) and window function (bottom) 
for ELODIE residual radial velocities of HD 196885 A when the 
contributions from the two bodies fitted in Tab|2]are subtracted. 
An important peak is detected at P — 368 days, which could 
be interpreted as a third body in the system. However, looking 
at the window function, we see that the same peak is present, 
indicating that this may be a false alarm. 

Curiously, performing a frequency analysis of the ELODIE 
radial velocity residuals (FigQ}, we find an important peak sig- 
nature at about P = 368 days (Fig|3]l, that could be interpreted as 
a second planet around HD 196885 A. However, this peak is not 
present when we analyze CORALIE data, casting some doubts on 
its origin. Computing the ELODIE window function, we also re- 
mark an important peak at exactly 368 days (Figf3]l, clearly sug- 
gesting that the peak shown in the frequency analysis is nothing 
but an aliasing of the observational data. Performing a Keplerian 
fit to the ELODIE residuals, we find that the best fit corresponds 
to solutions with very high eccentricities (e ~ 0.8), that are not 
dynamical stable. Moreover, when we plot a phase-folded of the 
residual radial velocities (Fig|4|i, we can see that our observa- 
tions only cover half of the orbit, preventing us from fully con- 
strain it. Finally, we performed a Keplerian fit of the three bod- 
ies also including data from CORALIE and CORAVEL. We found 
yffi = 1.460 and rms = 11.30m/s, which does not represent a 
substantial improvement with respect to the system with only 
two companions (Tabf2]i. The presence of a companion around 
368 days can thus be discarded. 

Besides the false alarm companion at 368 days, we may ask 
if there are other companions at different orbital periods. For that 
purpose we used a genetic algorithm, since we cannot clearly 
isolate any other peak in the frequency analysis of Figure [3] 
The inclusion of a third companion in the system allow us to 
reduce the V? U P t0 1-25 and the rms — 10.5 m/s. This rep- 
resents only a slightly better adjustment of the model to the 
data that can be justified as a natural consequence of increas- 
ing the number of free parameters. Moreover, identical adjust- 
ments can be obtained with many orbital periods, as different as 
1.3 days and 10 days, frequently with very high eccentricity val- 
ues. Therefore, no other companion can be conclusively detected 
in the residuals from the orbital solution listed in Table [2] 

5. Discussion and conclusion 

In this paper we report the presence of a planet orbiting the 
HD 196885 A star, with an orbital period of 1349 days. We also 
detect the presence of th e stellar companion HP 196 885 B, that 
was recently observed bv lChauvin et al . (2006, 2007). Its orbit is 
not completely determined, but we are able to provide some con- 
straints: P > 40 yr, a > 14 AU and M sin i > 0.28 M . It is possi- 
ble to find a large set of orbits with much longer orbital periods, 
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Fig. 4. Phase-folded ELODIE residual radial velocities for 
HD 196885 A when the contributions from the two bodies fitted 
in Tab|2] are subtracted. Data are superimposed on a Keplerian 
solution with e = and P = 368 days. The respective residuals 
as a function of Julian Date are displayed in the lower panel. We 
see that ELODIE data only cover half of the orbital phase, pre- 
venting us from completely constrain it. Moreover, the inclusion 
of this third body does not improve the fit significantly. We then 
conclude that the presence of a third body at P = 368 days is a 
false alarm. 



that fit the observational data as good as the minimum mass or- 
bits (Tab|3]l. In particular, we find a sol ution at P — 120 yr w ith 
M sin i = 0.6 M , the mass estimated bv lChauvin et al.1 (120071) . 

The HD 196885 A system was previously described as a star 
orbited by a planet at 386 days by the California & Carnegie 
Planet Search Team, detected with data acquired with the Lick 
survey. Although the system is reported on their webpage and 
many other places, no publication describing the system is yet 
known. The data acquired with the ELODIE survey also present 
some signal around 368 days, but a more detailed analysis shows 
that this signal corresponds to an aliasing of the observational 
data, probably due to the similarities with the Earth orbital pe- 
riod. We thus do not confirm the existence of the planet obtained 
with measurements form the Lick survey. 

The ELODIE radial velocity residuals after subtracting the 
signal from the two companions still presents a O - C (Fig[T|) 
slightly above the precision of the instrument. This suggests 
that the system may still hide additional planetary compan- 
ions. However, we were unable to find them, even when us- 
ing a genetic algorithm. The excess in the residuals can also 
be a contamination from the spectrum of the stellar compan- 
ion HD 196885 B. We have searched for its presence in our 
ELODIE data using multi-or der TODCOR, a two-dimensional 
cross-correlation algorithm dZucker et alj|2 004). but did not find 
anything convincing. We have also searched for correlations be- 
tween line bisectors and radial velocity or residuals. No signifi- 
cant correlation could be detected. These negative results do not 
allow us to formally discard the blend scenario, but they render 
it an unlikely possibility. 

The planet HD 196885 Ab is the fourth to be discov- 
ered in a close binary with a separation smaller than 25 AU. 
The other cases already known are Gliese86 dOueloz et al.l 
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20001: lLagrange et alJ 120061) . y Cep dHatzes et alj l2003t 



Neuhauser et all 120071) and HD41004A (Zuc ker et al.1 120041) . 



These systems are ideal for carrying out combined astrometric 
and radial-velocity observations to constrain the binary dynamic 
properties and the possible impact of a close binary companion 
on planet formation and evolution. 
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